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SNOW DEPTH AND SNOW EXTENT USING VHRR DATA FROM THE NOAA-2 SATELLITE 


David F. McGinnis, Jr., John A. Pritchard, and Donald R. Wiesnet 
National Environmental Satellite Service, NOAA, Washington, D.C. 


ABSTRACT. The NOAA-2 environmental satellite provides daily cov- 
erage of the Earth in the visible (0.6-0.7 um) and thermal (10.5- 
12.5 um) spectral bands. The ground resolution of the Very High 
Resolution Radiometer (VHRR) is 1 km at nadir. This improved re- 
solution in the visible permits more detailed observations of 
snow features than was possible with previous operational sat- 
ellites. A densitometer examination of a visible-band image from 
Feb. 11, 1973, which shows heavy snow cover in considerable de- 
tail over areas extending from Alabama to North Carolina, indi- 
cates that, in general, there is direct correlation between in- 
creasing brightness and increasing snow depths. Digitized re- 
flectance data from the study area were compared with prestorm 
bare-ground digitized reflectance data of Feb. 6, 1973, to deter- 
mine the relation of snow reflectivity to snow depths. A para- 
bolic regression analysis of greatest satellite brightness ver- 
sus greatest snow depth for 211 data pairs produced a correla- 
tion coefficient of 0.84. 


I. INTRODUCTION 


The brightness of snow received as reflected radiation by visible-band sensors aboard me- 
teorological satellites has been mentioned as a possible indicator of snow depth (Barnes and 
Bowley 1968, McClain 1973). Because of the much improved resolution of the Very High Reso- 
lution Radiometer (VHRR) data from the NOAA-2 environmental satellite, compared with the 
data available to previous investigators, a decision was made to test this hypothesis. Op- 
portunity for such a study occurred on Feb. 11, 1973, immediately following a heavy snow 
storm that struck the southeastern part of the United States. 


II. SOUTHEAST SNOW STORM OF FEBRUARY 9-10, 1973 


The winter of 1973 will be remembered well in the southeastern United States for the storm 
of February 9 and 10. This storm, the result of the interaction of two sharply contrasting 
air masses, produced record snow depths. An outbreak of very cold Arctic air entered the 
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Figure 1.--Total snowfall, snow storm of Feb. 1-10, 1973 


Northern Plains on Monday, Feb. 5, 1973. This large air mass plunged rapidly southward into 
Texas but made only slow progress eastward. By 1800 GMT (1300 LT) on Thursday, February 9, 

the leading edge of the cold air mass lay along the eastern slope of the Appalachian Mountains. 
Rain began falling in a broad band 350 km wide in the area of frontal convergence. Tempera- 
tures remained mild along the Atlantic seaboard until a developing cyclone drew in the very 
cold air just west of the Appalachians. The precipitation in Alabama and Georgia subsequently 
changed to freezing rain and sleet and then to snow during the morning and early afternoon of 
Friday, February 9. This snow spread into North Carolina and South Carolina as the deepening 
cyclonic disturbance moved northeastward into the Atlantic Ocean. Snow continued falling dur- 
ing most of Saturday in Georgia and the Carolinas, ending along most Carolina coastal sections 
after midnight. On Sunday morning, record snow depths were found on the ground at many loca- 
tions in Alabama, Georgia, and the Carolinas. Figure 1 shows the total snowfall in Alabama, 
Georgia, and South Carolina as analyzed by the individual State climatologists from scattered 
point observations taken largely by untrained cooperative observers. Data and isopleths, 
originally in inches, were obtained from individual State climatological summaries for Feb- 
ruary 1973 (Environmental Data Service 1973). Despite inconsistent values from State to State, 
particularly along the Georgia-South Carolina border (where the 2.5-cm isopleth in northeastern 
Georgia matches the 15-cm isopleth in northwestern South Carolina), values in excess of 30 cm 
were reported throughout an extensive band from east-central Alabama through northeastern 

South Carolina. The dashed line represents the authors' interpretation of the 2.5-cm isopleth 


in Alabama. Storm totals, because of the warm ground and settling of the snow, were somewhat 
greater than observations taken Sunday morning at 1200 GMT after the storm's ending. Rapid 


clearing followed the storm; and at approximately 1545 GMT (1045 LT) on orbit 1489, Feb. 11, 
1973, the NOAA-2 satellite viewed the cloud-free but snow-covered landscape (fig. 2). 


IIT. NOAA-2 SATELLITE 


The NOAA-2 satellite scans Earth in 25 orbits every 2 days at an altitude of approximately 
1500 km in a sun synchronous near-polar orbit. ! Essentially the same area of Earth is viewed 


"as of Dec. 17, 1974, NOAA 3 (with sensors identical to those on NOAA 2) was in a standby 
mode. NOAA 4 launched on Nov. 15, 1974, became the operational satellite on Dec. 17, 1974. 
a 


every other day, although the orbit of the sat- 
ellite progresses eastward about 1° of longi- 
tude every 25 orbits. Aboard the spacecraft 
are various sensors (Schwalb 1972), each with 
its own identical backup system. 


Most important of the sensors for hydrologi- 
cal use is the VHRR. This instrument scans 
Earth in both the visible (0.6-0.7 um) regions 
of the spectrum. The VHRR system generates 
analog signals that, after transmission to the 
readout station, are transformed into a photo- 
graphic format. VHRR data are also available 
on digital tapes. With a resolution of approx- 
imately 1 km at nadir, the VHRR imagery is 4 
times better in the visible and 10 times better 


in the thermal regions than data taken by in- 


Figure 2.--Portion of the NOAA-2 VHRR visible struments aboard NOAA 1 and earlier satellites. 
(0.5-0.6 um) image taken at 1545 GMT on 
Sunday, Feb. 11, 1973 
One of the drawbacks of VHRR images is that 


panoramic distortion increases with increasing distance from the center line of the image. 
The image is built up through a combination of the satellite's forward motion in orbit and 
its VHRR scanning perpendicular to this motion. This distortion is caused, in part, by the 
curvature of Earth's surface and produces foreshortening of the image toward the horizon. 


IV. IMAGERY ANALYSIS 


Figure 2 shows the snow-covered area as seen in the visible range by the VHRR. State bound- 
aries have been added to point out the distortion in the image and to provide reference points 
for location. These boundaries were added to figure 2 by using a Zoom Transfer Scope (ZTS), 
an optical device used to transfer information from unrectified satellite images to base maps, 
and vice versa. The distortion is greatest in the eastern portion of the image. A distinct 
western limit of the snow is evident, as well as a highly reflective (bright) strip in the 
center of the snow area, corresponding generally to the deepest snows shown in figure 1. Some 
land features in the snow belt, such as the swampy or heavily forested areas adjacent to many 
rivers, contrast sharply with adjacent areas in the snow-whitened land. The more heavily 
forested Piedmont Plateau in the western section of the snow area is mottled in appearance 
and is less bright than elsewhere. The 1-km resolution of the VHRR permits easy identifica- 
tion of many small rivers and also enables the analyst to determine the location of the west- 
ern snow edge to within 8 km when drainage features are used as local reference points. 


In figure 3 is a comparison of the snow extent viewed by the satellite with the 2.5-cm iso- 
pleth determined from ground observations on Sunday, February 11. Location of the zero iso- 
pleth, although the most useful for comparison with the satellite image, requires interpre- 
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Figure 3.--Comparison of snow extent as viewed by the NOAA-2 satellite with the 
2.5-cm isopleth determined from ground observations. The heavy solid line is 
the snow boundary based on the NOAA-2 VHRR visible imagery; and the dashed 
line is the snow boundary based on ground observations. 


tation of available ground truth. In areas where ground 


1.0 observations are few, accurate positioning of the zero 
isopleth is rather difficult. Fortunately, the sharp 
9 gradient in snow depths along the western boundary 
8 would place the edge of the snow only slightly north- 
west of the 2.5-cm isopleth. Note that Atlanta had no 
si 7 measurable snow during the storm although it lies less 
5 ¥ than 15 km from the 2.5-cm isopleth. Close agreement 
co exists between ground and satellite observations along 
= Ee) the western snow boundary. Even southeastern Georgia, 
a A where the edge of the snow line appears quite vague, a 
result of mixed precipitation and rather warm ground, 
3 satellite and ground maps of the snow extent coincide 
2 closely. 
| Analysis of the recorded brightnesses takes into ac- 
0 count the system's limitations. Each pixel (picture 
ne 6 af .8 element) of the VHRR image covers an area of approxi- 
WAVELENGTH (MICROMETERS } mately 1 kin (at nadir) on the ground. The brightness 
Figure 4.--NOAA-2 VHRR visible rel- value of each pixel thus represents an average response 
(eee anneal Borge ret RS that may include dark, forested regions as well as high- 


ly reflective open fields, even when both are snow covered. Also, the VHRR visible-band 
sensor is not a perfect detector. Figure 4 shows that the relative response of the detector 
falls short of perfect (1.0) throughout the design range of 0.6 to 0.7 um and also senses 
some energy outside these limits. 


The enhancement of land features by the differing reflectance of snow cover has been des- 
cribed. While this factor is useful for identifying and locating specific landmarks, it 
tends to mask the relation of snow brightness to snow depth. The aging of a snow pack is 
manifested in a reduction of the surface brightness--another factor to consider when attempt- 
ing to relate snow brightness to snow depth.. In our study of this relationship, the image 
was enhanced using a color densitometer. Various color schemes were tried, but none could 
discriminate more than three ranges of depth; these were 0 to 8 cm, 8 to 30 cm, and more 
than 30 cm. This breakdown into three ranges of depth is entirely qualitative and shows 
little improvement over the earlier estimates by Barnes and Bowley (1968) using lower reso- 
lution satellite imagery. One should note, however, that their study was confined to the 
Great Plains where variation in land use is much less of a problem. They were able to dis- 
tinguish only three distinct levels of brightness within snow-covered areas; these correspond 
to depths of a trace to 2.5 cm, 2.5 to 10 cm, and greater than 10 cm. 


In figure 5, a north-south densitometric trace through the snow area is shown on the left; 
the white strip through the image is the cross section that the trace represents. The densi- 
tometer trace is a quantitative measure of changes in photographic gray scales; it was pro- 
duced from a black and white photograph using a color densitometer. Brightness values are 
represented along the horizontal axis, with increasing values to the left. Outside the snow 
area, brightness values are uniformly low. Within the snow-covered area, values are higher 


but show several fluctuations, which are a function of varying land use. The abrupt termin- 
ation of the snow along the northwest boundary of the snow area is associated with a rapid 
decrease in brightness. The southern boundary of the snow area is more diffuse; thus, the 


decrease in brightness is gradual. 


V. DIGITAL TAPE ANALYSIS 


Digital computer printouts analogous to the pictorial display were produced from the VHRR 
analog tapes. The analog tapes were converted into nine-track digital tapes containing every 
scan line generated by the VHRR. Each scan line consists of individual, partially overlapping 
samples, two of which cover a 1-km? field of view. These elemental pieces of information can 
be averaged or combined in various schemes to meet the needs of the investigator. 


In digital printout, 36 different alphanumeric characters are available to cover the 255 
possible steps in the gray scale of the VHRR image. Usually, fewer than 255 steps are re- 
quired for analysis. In the case of the February 11 image, only 63 gray levels were used 


Figure 5.--Densitometric trace through portion of snow-covered area in South Carolina and 
Georgia 


within the snow area to represent brightness variations ranging from dark swamps to bright 
snow-covered cleared fields. For this study, a three line four-sample (approximately 3 km) 
moving average was selected. 


In pursuit of a snow depth-brightness relationship, seven brightness cross-sections of 
South Carolina, Georgia, and Alabama were plotted; each cross section was three scan lines 
(approximately 3 km) wide. A snow-free VHRR image (Feb. 5, 1973, orbit 1414) of this same 
area then was digitized to determine the brightness of the background with no snow. Then, 
the seven cross-sections from the snow-free data carefully were matched manually with those 
from the February 11 snow-depth data. 


The eastward precession of the NOAA-2 satellite orbit, about 1° longitude every 2 days, 
prevented an exact matching of the two orbits. Using landmarks as reference points, however, 
we achieved a close correspondence. Figure 6 shows one of the seven matched pairs of cross 
sections. Note that the snow-free terrain is nearly the same brightness everywhere. The 
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zero values indicate the location of one of the two subsynchronization lines that are part 
of the VHRR data stream. The same swath under a blanket of snow varies widely in brightness, 
demonstrating how surface features are brought out by the snow. Subtraction of the snow-free 
terrain brightness from the snow brightness did not provide improved insight into the snow 


brightness-snow depth relationship. 


Qwing to the scarcity of detailed snow-depth measurements on the morning of Feb. 11, 1973, 
comparison of the snow brightness with ground truth necessitated substitution of snowfall 
totals [as provided in Environmental Data Service (1973)] for snow depth values. Arbitrary 
values between 1 and 256, corresponding to the 256 brightness increments detectable by the 
VHRR, were assigned to the 3-km data points in the seven cross-sections. A simple linear 
regression relating these brightness values to snowfall totals (assumed to vary linearly be- 
tween the isopleths) was run for each cross section in another attempt to determine how well 
brightness relates to snow depth. The correlation coefficients ranged from 0.10 to 0.20, 
disappointingly low. These values probably were due to the use of 3-km averaged data, in 
which varying types of land use, some brighter than others,are combined to produce an aver- 
aged brightness value. Varying land use, when enhanced by snow, provides land marks for 1lo- 


cation but masks the snow-brightness to snow-depth relationship. 
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Figure 6.--Visible reflectance values before and after the snow 
storm for a 3-km cross section in Georgia and Alabama 
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A computer program was adapted to accept VHRR digital data and produce, at 1-km or 4-km 
resolution, histograms of brightness values at preset class intervals and provide the small- 
est and greatest brightness values for each 32 x 32 element area. The program also contours 


each area using the preset class intervals, thus providing an isopleth map. The 1,024 bright- 
ness values and annotations for location are included in the program. 


Using the 1-km resolution option (32 x 32 km squares), we fitted a grid of 201 VHRR squares 
to figure 1. The greatest snow depth for each square (from fig. 1) was matched with the 
greatest brightness value. Linear and parabolic best fits were calculated for the data pairs 
with brightness as the dependent variable. Figure 7 shows the 201 data points and the best 
fitting linear and parabolic curves. Correlation coefficients were 0.80 for the straight 
line and 0.84 for the parabolic curve. 


The parabolic fit of the data appears more consistent with the apparent relationship of 
brightness to snow depth. Rather small increases in depth when the snow is only several cen- 
timeters deep produce large increases in brightness. In this case, the snow covers more and 
more objects on the ground as it deepens, first individual blades of dormant grass, then 
small shrubs, then larger objects. Once the snow accumulates to about 25 cm, most small 
plants are covered, only the larger shrubs and plants (corn stalks to trees) remain visible 
and in most cases will remain uncovered except in some mountainous areas where extremely 
heavy snows occur. 


Some radiation penetrating into a snow pack, if not absorbed by underlying vegetation and 
soil, eventually returns to the surface and contributes to the albedo. Giddings and LaChapelle 
(1961) demonstrated that absorption from an underlying black surface greatly reduced the al- 
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Figure 7.--Relation of relative VHRR visible reflectance to depth of snow 
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bedo of aged snow until snow depths exceeded 20 cm. Thus small increases in satellite- 
observed brightness may be observed even when tops of vegetation are being covered with 10 
to 20 cm of snow. 


Regarding data discontinuities in the snow isopleths at State boundaries, we used the fol- 
lowing procedure. When a 32-km square fell across a State boundary, the greatest depth was 
determined using data from the State having the greater portion of the square. Also, some 
error in snow-depth estimates could occur because of the misalinement of distortion of the 
Squares fitted to figure 1. Further, most of the snow-covered area in Alabama received some 
sunshine late in the day on Saturday, 10 February 1973, which probably melted a small amount 
of snow. The amount of snow remaining at the time of the satellite observation on Sunday 
morning, therefore, is believed to be slightly less than that indicated in figure 1. The 
scarcity of ground observations precluded a re-evaluation of snow depths. The effect of the 
Sun on the reflectivity of the snow and, thus, on the brightness recorded by the satellite 
also is not considered. Both of these effects caused by the Sun, however, trend in the same 
direction (i.e., toward decreasing values). 


The curves derived from the data are for the NOAA-2 VHRR visible-sensor 1. Equivalent 
curves may be different for either NOAA-3 or NOAA-4 visible sensors or the other NOAA-2 VHRR 
visible sensor. The general relationship likely would remain, but with perhaps differing 
brightness limits. For this case, the lower limit of brightness (bare ground) is between 40 
and 43; the upper limit of snow brightness is 101. 


VI. RESULTS 


The improved resolution of the VHRR over previous sensors on operational satellites permits 
mapping of snow limits to better than a 10-km resolution, which is an improvement over the 
20- to 40-km resolution reported with earlier satellites (McClain 1973). Although the radi- 
ance data received from the satellite are unrectified, the enhancement of identifiable land 
features by the snow permits the precise location of snow boundaries of mesoscale resolution. 
Color renditions of the black and white imagery by means of a color densitometer, although 
subjective, do permit coarse estimates of snow depths, but the improvement over estimates 
made from earlier, lower resolution satellite sensors is minor. 


Linear regression analyses of cross sections relating snow brightness to storm totals pro- 
duced negligible correlation coefficients. A major reason for the lack of correlation is the 
variable brightness of snow-covered land used for different purposes. For snow-free winter 
conditions, all types of land surfaces in the southeastern United States appear almost uni- 
formly bright (lower curve in fig. 6). When enhanced by a deep snow cover, however, the vary- 
ing land surfaces reflect very differently (upper curve in fig. 6). Under the same depth of 
snow, open farmland appears very bright, wet swampy areas adjacent to the major rivers are 
dark, and forested areas appear in differing intermediate tones depending upon density and 


type of trees. 


A significant correlation of 0.84 was achieved, however, when the greatest brightness value 
in a given 32 x 32 km square was paired with the greatest snow depth for the same square. 
Using the greatest brightness within a particular region is, in effect, searching for the 
clearings and open areas that are free of trees and other nonsnow-covered objects that lower 
the average brightness of the 1-km element. 


VII. CONCLUSIONS 


For the case studied (i.e., freshly fallen snow over an area of widely varied land use), 
the VHRR provides a synoptic picture of snow extent not possible from point observations on 
the ground. The 1-km resolution of the VHRR sensors permits positioning of snow limits to 
within 10 km, some two to four times better than that possible with data from earlier opera- 


tional satellites. 


Snow-depth estimates correlate well with snow brightness. The parabolic relationship indi- 
cates that a rapid increase in brightness corresponds to increases in snow depths when depths 
are less than 25 cm; for depths exceeding 25 cm, much smaller brightness increases are en- 


countered as depths increase. 
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